First Results of the Phase II SIMPLE Dark Matter Search 
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We report results of a 14.1 kgd measurement with 15 superheated droplet detectors of total active 
mass 0.208 kg, comprising the first stage of a 30 kgd Phase II experiment. In combination with the 
results of the neutron-spin sensitive XENON10 experiment, these results yield a limit of \a p \ < 0.32 
for Mw = 50 GeV/c 2 on the spin-dependent sector of weakly interacting massive particle- nucleus 
interactions with a 50% reduction in the previously allowed region of the phase space formerly 
defined by XENON, KIMS and PICASSO. In the spin-independent sector, a limit of 2.3xl0~ 5 pb 
at M w = 45 GeV/c 2 is obtained. 



The direct search for weakly interacting massive par- 
ticle (WIMP) dark matter continues to be among the 
forefront endeavors of modern physics activity. Search 
experiments are generically based on the detection of nu- 
clear recoil events resulting from WIMP-nucleus interac- 
tions, and are traditionally classified as spin-independent 
(SI) or spin-dependent (SD) according to which interac- 
tion channel the experiment is most sensitive, of which 
the first has generally attracted the most attention. The 
current status of the SI search for WIMPs is defined by a 
number ofprojects, including XENON Q, CDMS and 
ZEPLIN [3j, which as a result of their target nuclei spins 
also define the WIMP-neutron sector of the SD phase 
space. The WIMP-proton sector is currently constrained 
by PICASSO [4] and KIMS @. 

The SIMPLE (Superheated Instrument for Massive 
ParticLe Experiments) [1,01 project, located in a 61 m 3 
cavern at the 1500 mwe level of the Laboratoire Souter- 
rain a Bas Bruit (LSBB) in southern France, currently 
runs superheated droplet detectors (SDDs). The SDD is 
a suspension of 1-2% superheated liquid C2CIF5 droplets 
(~ 30 /zm radius) in a viscoelastic 900 ml gel matrix 
which undergo transitions to the gas phase upon en- 
ergy deposition by incident radiation. Two conditions 
are required for the nucleation of the gas phase in the 
superheated liquid @: (i) the energy deposited must 
be greater than a thermodynamically-dcfined minimum 
energy, and (ii) this energy must be deposited within 
a thermodynamically-defined maximum distance inside 
the droplet. Together, energy depositions of order ~ 
150 keV/^m are required for a bubble nucleation, which 
renders the SDD effectively insensitive to the major- 
ity of traditional detector backgrounds which complicate 
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more conventional dark matter search detectors (includ- 
ing electrons, 7's and cosmic muons). The insensitivity 
is not trivial, comprising an intrinsic rejection factor su- 
perior to that of other search techniques by 1-5 orders 
of magnitude. Additional advantages of the superheated 
technique include low cost, scalability, and increased sen- 
sitivity to WIMP-proton spin interactions via the 19 F 
content [Toj . 

The SDDs were fabricated according to previously- 
described procedures @, in an underground (210 mwe) 
clean room in close proximity to the measurement 
site. The SIMPLE gel ingredients, all biologically-clean 
food products, are purified using actinide-specific ion- 
exchanging resins. The freon is single distilled; the water, 
double distilled. The presence of U/Th contaminations 
in the gel, measured at ~ 0.1 ppb by low-level a and 7 
spectroscopy of the production gel, yields an overall a- 
background level of < 0.5 evt/kg freon/d. A similar level 
is measured for the detector containment materials. 

The detectors are capped using a mechanical construc- 
tion which virtually eliminates pressure microleaks Q- 
Each cap contains feedthroughs for pressure monitor- 
ing, and for a 20-16k Hz electret microphone encased 
in a latex sheath, which is immersed in a 4 cm thick 
glycerin layer covering the gel at the top of the detec- 
tor containment. Each microphone response is pream- 
plified and recorded in a MatLab platform in sequen- 
tial files of 8 min duration [H, with resolutions of 0.3 
mV in amplitude and 1.6 x 10~ 2 ms in time; the pres- 
sure reading is similarly recorded separately. The use 
of shielded telecommunications-grade cabling eliminates 
signal resulting from cable motion and parasitics, even 
when exaggerated. 

The SDDs are immersed to a depth of 20 cm in a 
700 liter water pool maintained at a bath temperature of 
9.0±0.1 °C within the cavern, and pressurized to 2 bar to 
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reduce background sensitivity. The water pool rests on a 
dual vibration absorber placed atop a 20 cm thick wood 
platform resting on a 50 cm thick concrete floor. The 
pool is surrounded by three layers of sound and thermal 
insulation. An additional 50-75 cm thick water shielding 
surrounds the insulated pool and pedestal, with a 75 cm 
water thickness overhead; 50 cm of water separates the 
pool bottom from the SDD bases. 

At 1500 mwe, the ambient neutron flux is primarily due 
to the surrounding calcite rock, estimated at well-below 
4 x 10^ 5 n/cm 2 s The cavern is shielded from the 

rock environment by a 30-100 cm thickness of concrete, 
internally sheathed by a 1 cm thickness of iron. Radio- 
assays of the concrete yielded 1.90±0.05 ppm 232 Th and 
0.850±0.081 ppm 238 U; of the steel, 3.20±0.25 ppb 232 Th 
and 2.9±0.2 ppb 238 U. The results are at the same level 
as those recorded in other underground locations such as 
Canfranc, Modane and Gran Sasso [lH-Q- Monte-Carlo 
simulations of the on-detector neutron field, which in- 
clude all shielding materials and account for spontaneous 
fission plus decay-induced (a,n) reactions, show negligi- 
ble variations for concrete thicknesses > 20 cm, and yield 
an expected neutron background of 1.09 ± 0.02 (stat) ± 
0.07 (syst) evt/kgd. 

The ambient radon level varies seasonally between 28 - 
1000 Bq/m 3 as a result of water circulation in the moun- 
tain. The cavern air is purged ~ 10 times per day, re- 
ducing the ambient radon levels to ~ 60 Bq/m 3 . Diffu- 
sion of the environmental radon into a detector is limited 
by the surrounding water, which covered the detectors 
to just above their glycerin levels, and is circulated at 
25 liter/min (equivalent to replacing the top 1 cm wa- 
ter layer each minute). The radon contribution is also 
low because of the short radon diffusion lengths of the 
SDD construction materials (glass, plastic, metal), the 
N2 overpressuring which inhibits the advective influx of 
Rn (via stiffening of the gel), and the glycerin layer cov- 
ering the gel. The reduction in the overall radon con- 
tribution to the measurement, including its progeny, is 
estimated at ~ 10 5 with the overall a contribution to 
the measurement (including the detector contribution) 
estimated at 3.26 ± 0.08 (stat) ± 0.76 (syst) evt/kgd. 

Data obtained from 15 SDDs, containing between 8-21 
g of C2CIF5 for a total active mass of 0.208 kg, between 
27 October 2009 and 05 February 2010 were analyzed for 
this report. An additional, similarly installed, freon-less 
but otherwise identical SDD, served as an acoustic veto. 
The total raw exposure was 14.10±0.01 kgd, with 1.94 
kgd resulting from the detectors being introduced at one 
device per day over the three week installation period, 
and a 4.70 kgd loss from weather-induced power failures 
during the run. 

The SDD signals, pressures and temperature are mon- 
itored continuously during operation, as also the radon 
level. Each detector was first inspected for raw signal rate 
and pressure evolution over the measurement period. An 
initial data set (4056 events) was then formed by pass- 
ing the data files through a pulse validation routine [8[ 



which tagged signal events if their amplitude exceeded 
the noise level of the detector by 2 mV. Tagged signals 
in coincidence with the freon-less device were next re- 
jected, as also all candidate signals with less than five 
pulse spikes above threshold; the remaining set was then 
cross-correlated in time between all SDDs, and coinci- 
dences rejected as local noise events and that a WIMP 
interacts with no more than one of the in-bath detectors. 
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FIG. 1: Scatter plot of the squared amplitude and frequency 
of the primary harmonic of each true nucleation event (solid), 
together with the same for neutron and a calibration events 
(open), with the dashed horizontal line corresponding to a 
signal amplitude of 100 mV. 

The signal waveform, decay time constant and spectral 
density structure of the remaining 1828 single events were 
next inspected individually. A particle-induced nucle- 
ation event possesses a characteristic frequency response, 
with a time span of a few milliseconds, a decay constant 
of 5-40 ms, and a primary harmonic between 0.45-0.75 
kHz This response differs significantly from those of 
gel-associated acoustic backgrounds such as trapped N2 
gas (3.4% of the total), gas escape (0.008%), and gel frac- 
tures (4.4%) which appear at lower frequencies [8|, as well 
as local acoustic backgrounds (88%) such as water bub- 
bles which differ in power spectra. This event-by-event 
analysis permits isolation of the true nucleation events 
with an efficiency of better than 97% at 95% C.L. Fig. 
1 displays the signal amplitudes (A) with frequency for 
each of the identified 60 particle-induced signal events, 
in which a gap corresponding to A = 100-130 mV is dis- 
cernible. 

At 9°C, the reduced superheat of the devices is 0.3, 
and the probability of events from electrons, 7's and 
mip's neglig ible over the exposure. Calibrations of the 
a response have been made by doping the devices with 
U308 during fabrication. The event signals, identified in 
the same fashion as described above, are shown in Fig. 
1. The SRIM-calculated dE/dx for a's in C 2 C1F 5 has a 
Bragg peak which sets a lower energy threshold at 200 
keV for the operating temperature and pressure; below 
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this threshold, ev's are detected only through a-induced 
nuclear recoils. 

Calibration of the high concentration SDD response to 
neutrons, using sources of Am/Be, yielded a minimum 
threshold recoil energy (E t / lr ) of 8.0±0.1 keV, with an 
acoustic detection efficiency of 0.98±0.03 at 9°C and 2 
bar. The events are also displayed in Fig. 1, all of which 
occur with amplitudes < 100 mV. The difference in the 
two distributions, particularly at lower amplitudes, re- 
sults from performing the calibrations with a 15 cm wa- 
ter shield to enhance the tails on the moderated neutron 
spectrum. 

Fig. 2 displays a typical histogram of both calibration 
signal amplitudes. As seen, the neutron-induced events 
are of lower amplitude than the a-induced, and empir- 
ically well fit by a Gaussian plus constant background 
from which a discrimination cut for A < 100 mV is placed 
with an acceptance of > 97%. The small droplet size pro- 
vides a natural lower cutoff to the deposited a energy and 
is together with the dE/dx responsible both for the am- 
plitude gap between the a and neutron populations, and 
the spectral asymmetry in the a distribution. 
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FIG. 3: a p -a„ for SIMPLE at M w = 50 GeV/c 2 , together 
with benchmark experiment results; the dashed contour rep- 
resents a "0 evt" SIMPLE result for the same exposure. The 
allowed regions are defined by a band (single nuclei target) or 
ellipse (multinuclei target), with the external area excluded. 
The cross-hatched central area about (0,0) indicates the re- 
gion allowed by this result and XENON. 

The impact of the result in the SD phase space at M\y 
= 50 GeV/c 2 is shown in Fig. 3. The contour is cal- 
culated within a mo del- independent formulation (Iol |. in 
which the region excluded by an experiment lies outside 
the indicated band, and the allowed region is defined by 
the intersection of the various bands. The calculations 
use a standard isothermal halo, bubble nucleation effi- 
ciency of (1-E thr /E) 0, and WIMP scattering rate [H 
with zero momentum transfer spin-dependent cross sec- 
tion ijsd for elastic scattering: 



FIG. 2: A typical histogram of calibration neutron and a am- 
plitudes, each population consisting of several hundred events. 
The vertical dashed line indicates the neutron discrimination 
cut at A < 100 mV, which includes all neutron calibration 
events. 

The 14 low amplitude events of the run are consistent 
with the neutron calibrations, yielding 0.99 ± 0.27 (stat) 
evt/kgd. Corrected for detection and identification ef- 
ficiencies, the difference between measurement and neu- 
tron background estimate is better than 0.3<t. 

An upper limit on the number of WIMP events in 
the presence of the uncertain neutron background is es- 
timated by applying the Feldman- Cousins method [llij . 
based on observing 14 events with a background one stan- 
dard deviation below the central value of the expected 
neutron background. This yields 4.3 events, and a re- 
sulting WIMP rate of 0.57 evt/kgd at 90% C.L. 



(j S d ~ G F [a p < S p > +a n < S n >] 2 - , (1) 

where et p .„ are the WIMP-proton, neutron coupling 
strengths, < S p . n > are the expectation values of the 
proton (neutron) group's spin, Gf is the Fermi coupling 
constant, and J is the total nuclear spin. The form fac- 
tors of Ref. have been used for all odd-A nuclei. The 
spin values of Strottman have been used for 19 F [l8|; for 
35 C1 and 37 C1, < S p , n > are from Ref. 0, while for 13 C 
the < S Pt n > were estimated by using the odd group ap- 
proximation. Use of the Divari et. al. spin values [19j for 
19 F would rotate the ellipse about the origin to a more 
horizontal position. We include only the 121 kgd result of 
CDMS since the more recent [2(| mo del- independent 
result hasn't yet been published. The shaded area rep- 
resents the allowed DAMA/LIBRA region [2l|, which is 
already excluded by other experiments. 
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As indicated, the present result combined with 
XENON10 yields limits of \a p \ < 0.32, |a„ < 0.16 on 
the SD sector of WIMP-nucleus interactions for 
GeV/c 2 , with ~ 50% reduction in the allowed region of 
the phase space. M\y above or below this choice yield 
less restrictive limits 11011 . 
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FIG. 4: Spin-independent contour for SIMPLE; the dashed 
contour represents a "0 evt" SIMPLE result for the same mea- 
surement exposure. Also shown are several [2^] of the leading 
spin-independent search results; the shaded area represents 
the recent result of CoGeNT [23[. 

For completeness, the impact of the result in the SI 
sector, calculated following the standard isothermal halo 
and WIMP elastic scattering rate of Ref. [T^J using a 
Helm form factor, is shown in Fig. 4 in comparison with 
results from other leading search efforts [22[. The fig- 
ure indicates a contour minimum of 2.3 xlO -5 pb at 45 



GeV/c 2 , slightly less restrictive than that of the 52 kgd 
COUPP [24[ exposure. The near-equivalence of the two, 
despite the ~ 48 x difference in sensitivity from coherence 
enhancement, most likely derives from the high radiopu- 
rity and a-discrimination of the SIMPLE detectors. 

In summary, a conservative analysis of 14.1 kgd of 
data from the first phase run of SIMPLE yields new re- 
strictions on the WIMP-proton parameter space of SD 
WIMP interactions which improves on those from both 
PICASSO and COUPP. These results represent the first 
stage of the ongoing 30 kgd Phase II exposure permit- 
ted by current funding, and further demonstrates the 
competitiveness of the superheated liquid technique in 
the search for astroparticle dark matter in both spin 
-dependent and -independent sectors. To further ap- 
proach the dashed contours of SIMPLE in Figs. 3 and 
4 requires further reduction and/or elimination of the 
neutron background: improvements in the second stage 
of this measurement include a 10 cm increase in the 
wood pedestal beneath the waterpool, and installation 
of an additional 10 cm polyethylene beneath the detec- 
tors, which together are projected to reduce the overall 
neutron contribution by ~ factor 5. This is accompanied 
by additional a and neutron calibrations, and a refined a- 
neutron discrimination analysis (25j | towards understand- 
ing the differences between these results and those of PI- 
CASSO nu. 
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